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ABSTRACT 

A one-dimensional  model  has  been  used  to  investi- 
gate the  effect  of  inaccuracies  in  meteorological  para- 
meters on  forecasts  of  upper  ocean  thermal  structure.  Ver- 
tical eddy  fluxes  of  heat,  salinity,  and  momentum  in 
the  oceanic  mixed  layer  are  parameterized  using  the  Level- 
2 turbulence  closure  model  of  Mel  lor  and  Yatnada  (1974). 
Surface  eddy  fluxes  of  sensible  heat,  latent  heat,  and 
momentum  are  related  to  the  imposed  meteorological  para- 
meters through  a bulk  aerodynamic  formulation.  Surface 
fluxes  of  solar  and  infrared  radiation  are  paramaterized 
with  the  formulas  of  Wyrtki  (1965). 

The  temperature  anomaly  in  the  ocean  forecast  pro- 
duced by  inaccuracies  in  the  meteorological  parameters 
diffuses  rapidly  downward  from  the  surface  to  within  a few 
meters  of  the  mixed  layer  base.  Its  downward  diffusion 
is  much  slower  below  this  level  and  this  tends  to  produce 
a relatively  large  vertical  gradient  in  temperature  anomaly 
near  the  base  of  the  mixed  layer  which  may  have  important 
implications  for  acoustic  modeling.  Both  the  depth  of  the 
mixed  layer  and  the  magnitude  of  the  mixed  layer  temperature 
anomaly  are  found  to  be  nonlinear  functions  of  the  in- 
accuracies in  the  potential  temperature,  humidity,  and  wind 
speed  of  the  atmospheric  mixed  layer. 


Estimates  of  the  time  evolution  of  the  inaccuracies 
in  the  air  temperature,  humidity,  and  wind  speed  during 
an  actual  operational  forecast  situation  are  made.  For 
the  typical  summertime  subtropical  stratification  considered 
here,  the  temperature  anomalies  resulting  at  the  end  of 
a five  day  forecast  from  these  inaccuracies  alone  are  a 
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INTRODUCTION 


The  sophistication  of  the  Navy's  various  acoustic 
detection  systems  has  increased  in  recent  years  and  made 
detailed  knowledge  of  the  oceanic  environment  of  prime 
importance.  Presently,  however,  the  Navy  his  no  capability 
to  forecast  the  internal  structure  of  the  ocean  and  must  rely 
on  coarse  analyses  of  historic  and  synoptic  data  for  its 
needs.  The  Naval  Oceanographic  Prediction  System  (OPS)  has 
been  proposed  to  provide  the  Navy  with  a better  knowledge  of 
the  oceanic  environment  on  an  operational  basis.  The  goal  of 
OPS  is  to  provide  five-day  predictions  of  the  major  features 
of  the  world  ocean  with  reasonable  accuracy  and  detail  and 
two-day  predictions  for  smaller  regions  with  higher  accuracy 
and  detail.  When  it  is  completed,  OPS  will  become  part  of 
the  Primary  Environmental  Processing  System  (PEPS)  at  Fleet 
Numerical  Weather  Central.  There  it  will  interface  with  a 
data  analysis  and  initialization  system  and  an  atmospheric 
prediction  model  to  provide  a dynamic  description  of  the 
marine  environment. 

Computed  ocean  forecasts  will  be  subject  to  error. 
The  reasons  for  these  forecast  errors  fall  into  three  basic 
categories  as  described  below. 

1)  The  instantaneous  state  of  the  ocean  is  never 
known  exactly.  Thus,  the  initial  conditions 
supplied  to  the  ocean  forecast  model  will  differ 
from  the  true  initial  conditions. 

2)  The  ocean  forecast  model  cannot  reproduce  the 
internal  dynamics  of  the  ocean  exactly.  This  is 


caused  primarily  by  the  fact  tu«*t  the  model  is 
discrete  in  time  and  space  and  cannot  resolve  all 
scales  of  motion. 

3)  The  state  of  the  atmosphere  and  therefore  the 
meteorological  forcing  of  the  ocean  are  not  known 
exactly  during  the  forecast  period.  This  is 
because  the  forecast  state  of  the  atmosphere  is 
also  subject  to  inaccuracies,  mainly  resulting  from 
inexact  knowledge  of  initial  conditions  and  incor- 
rect representation  of  internal  atmospheric  dynamics. 
An  additional  source  of  error  is  that  the  forecast 
sea  surface  temperature,  which  provides  important 
lower  boundary  conditions  for  the  model  atmosphere, 
will  not  be  perfect. 

In  this  paper  we  will  investigate  ocean  forecast 
errors  due  to  the  reasons  listed  in  category  3 above  by  inves- 
tigating the  sensitivity  of  5-day  forecasts  to  inaccuracies 
in  cloud  cover,  temperature,  humidity,  and  wind  speed  above 
the  sea  surface.  We  will  restrict  our  attention  to  regions 
of  relatively  high  horizontal  uniformity  and  will  concern  our- 
selves primarily  with  forecasting  the  vertical  thermodynamic 
structure  of  the  upper  ocean.  Thus,  this  study  can  be 
accomplished  with  a one-dimensional  model. 

Almost  universally,  the  upper  ocean  is  characterized 
by  a mixed  layer  extending  from  the  surface  to  about  5-100  m 
depth  in  which  temperature,  salinity,  and  current  velocity 
exhibit  only  a very  small  change  with  depth.  This  layer  owes 
its  high  degree  of  vertical  homogeneity  to  mixing  caused  by 
turbulence  that  is  generated  by  the  wind  and  intermittent 


upward  buoyancy  flux  through  the  sea  surface.  A dynamically 
stable  water  mass  in  which  the  vertical  eddy  fluxes  are  very 
small  exists  below  the  mixed  layer.  During  periods  of 
relatively  strong  wind  forcing  and/or  strong  surface  cooling, 
the  mixed  layer  tends  to  deepen  because  the  more  dense  water 
from  below  is  entrained  into  the  layer  by  turbulence  at  its 
base.  During  periods  of  relatively  weak  wind  forcing  and/or 
strong  surface  heating,  however,  the  source  of  turbulent 
kinetic  energy  may  become  too  weak  to  maintain  active  entrain- 
ment at  the  mixed  layer  base,  causing  the  layer  to  retreat  to 
a shallower  depth.  Thus,  the  problem  of  modeling  the  upper 
ocean  is  closely  associated  with  parameterization  of  turbulent 
processes  in  the  mixed  layer. 

The  first  attempts  at  modeling  the  upper  ocean  were 
with  diffusive  models.  In  this  type  of  approach,  vertical 
turbulent  fluxes  are  parameterized  in  terms  of  the  mean  fields 
using  eddy  diffusion  coefficients  which  may  themselves  depend 
on  the  mean  fields.  An  advantage  of  this  type  of  model  is 
that  no  assumptions  (other  than  the  eddy  coefficient  assump- 
tion) need  to  be  made  concerning  the  behavior  of  the  mixed 
layer . 


The  first  diffusive  model  that  predicted  the  density 
and  current  structure  of  the  upper  ocean  was  that  of  Munk  and 
Anderson  (1948).  The  forms  of  the  eddy  coefficients  in  this 
model  were  determined  empirically  from  stratified  turbulence 
data.  Mamayev  (1958)  and  Pandolfo  (1969)  have  also  proposed 
empirical  forms  for  the  eddy  diffusion  coefficients.  More 
recently,  models  such  as  those  of  Vager  and  Zilit inkevich 
(1968)  and  Mellor  and  Yamada  (1974)  have  been  developed  in 
light  of  modern  turbulence  closure  theory.  Martin  (1976)  and 
Jacobs  (1978)  have  performed  numerical  experiments  comparing 
several  of  the  proposed  forms  of  the  eddy  coefficients. 
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We  will  use  the  diffusive  (i.e.,  Level  2)  model  of 
Mellor  and  Yamada  (1974)  in  this  study.  This  model  has  been 
applied  to  the  upper  ocean  with  good  success  by  Mellor  and 
Durbin  (1975).  Martin  (1976)  concluded  that  it  was  the  most 
suitable  model  for  simulating  the  upper  ocean  among  those  he 
compared  and  it  has  been  proposed  as  an  important  component  of 
OPS  (Grabowski  and  Roberts,  1978). 

We  will  parameterize  the  surface  eddy  fluxes  of 
latent  and  sensible  heat  using  the  bulk  aerodynamic  method  in 
which  the  fluxes  are  related  to  mean  field  quantities  by 
means  of  bulk  exchange  coefficients.  This  technique  is  based 
primarly  on  the  theory  of  Monin  and  Obukhov  (1954)  and  has 
been  reviewed  comprehensively  by  Businger  (1975)  and  Busch 
(1977). 

The  remainder  of  this  paper  is  organized  as  follows: 
Section  2 describes  the  model  ocean.  Section  3 discusses  the 
parameterization  of  the  surface  fluxes.  Section  4 describes 
the  design  of  the  experiments.  Section  5 presents  the 
results.  Section  6 is  a summary. 
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2.0  DESCRIPTION  OF  THE  MODEL  OCEAN 

2 . 1 Fundamental  Equations 

1 

Assuming  zero  mean  vertical  motion  and  horizontal 
uniformity  for  all  fields,  the  conservation  equations  for 
heat,  salt,  and  momentum  in  the  upper  ocean  can  be  written 


9T  _9  /r  9'f\ 
9t  9z  y 11  dz ) 


If  = ^ (Ks  || ) 


^ X *•  _ g ( v 

St  - f * x I + si 


where  T is  the  temperature,  S the  salinity,  v the  horizontal 
current  velocity  vector  (standard  right-handed  Cartessian 
coordinate  system  with  z positive  upward  and  origin  at  sea 
surface  assumed),  F the  downward  flux  of  solar  radiation,  and 
f the  Coriolis  parameter.  Other  symbols  are  defined  in  the 
Appendix  and  all  the  notation  is  standard.  Spatial  averages 
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at  constant  depth  taken  across  a region  large  enough  to  encom- 
pass a statistically  significant  sampling  of  all  unresolved 
phenomena  are  denoted  by  (~)  and  implicit  in  (l)-(3)  are  the 
eddy  coefficient  assumptions 


-KhS 


-Kq  as 

s a z 


w'y ' 
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where  w is  the  vertical  component  of  the  motion,  and  the 
primes  indicate  departures  from  the  horizontal  averages.  The 
quantities  Kjj , Ks»  and  Ky  are  the  eddy  diffusion  coefficients 
and  are  determined  from  the  considerations  in  the  next 
section.  Finally,  we  note  that  all  horizontal  eddy  fluxes 
are  assumed  zero. 

2.2  Determination  of  the  Eddy  Diffusion  Coefficients 


Through  systematic  scaling  of  the  equations  for  the 
Reynolds  fluxes  and  turbulent  kinetic  energy,  Mellor  and 
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Yamada  (1974)  developed  a set  of  turbulence  closure  models 
for  planetary  boundary  layers  labeled,  in  order  of  increasing 
complexity,  Levels  1 through  4.  They  used  hypotheses  proposed 
by  Kolmogoroff  (1942)  and  Rotta  (1951)  to  model  the  triple 
correlation  and  dissipation  terms  in  these  equations  and  the 
empirical  constants  arising  from  the  theory  were  determined 
from  neutrally  stratified  turbulence  data. 


In  the  Level  2 model,  which  will  be  used  here, 
(4)-(6)  are  recovered  with  the  eddy  diffusion  coefficients 
given  by 


KH  = l(2Ef  Ah 


(7) 


Km  » 1(2E)*  AM 


(8) 


where  1 is  a turbulence  length  scale  and  E is  the  turbulent 
kinetic  energy 


E = |(u'2  + v'2  + w'2) 


(9) 


The  quantities  A^  and  Ay  are  stability  functions  given  by 
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where  the  n^-ng  are  empirical  constants  and 
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Note  here  that  the  mean  field  density  p in  (12)  can  be  related 
to  the  mean  field  temperature  and  salinity  with  a linearized 
equation  of  state 


p — p v/  1— a(T— T^/)  + b(S— Sy) 


(13) 


where  Pw  is  a reference  density,  Tw  a reference  temperature, 

Sw  a reference  salinity  and  a and  b empirical  constants.  Note 
also  that,  following  many  others,  we  assume 


Ks  = KH 


(14) 


The  turbulent  kinetic  energy  E is  calculated  from 
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The  three  terms  on  the  right  of  (15)  represent  shear 
production,  buoyant  damping,  and  viscous  dissipation  of  tur- 
bulent kinetic  energy  respectively. 

Finally,  an  estimate  of  the  turbulence  length  scale 
1 is  calculated  from  the  vertical  extent  of  the  turbulence 
field  according  to 


1 = 


r -» 

n8  J- « zE 


dz 


/’■* 

— CO 


(16) 


dz 


which  closes  the  model. 


Note  that  the  "mixed  layer"  is  henceforth  defined  as 
the  region  where  E>0.  Below  the  mixed  layer  we  retain  weak 
vertical  diffusion  by  setting 


Kg  = Kg  = Kf,j  = 1 cm^  S-1 


(17) 


For  a detailed  discussion  of  the  scaling  which  leads  to  this 
turbulence  closure  theory,  see  Mellor  and  Yamada  (1974). 
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PARAMETERIZATION  OF  THE  SURFACE  FLUXES 


The  total  fluxes  of  heat,  salinity,  and  momentum  at 
the  ocean  surface  are  given  by 


pwc(WT')0  = Fc  + B0  + H0  + LvQ0  (18) 


Pw(w'S')0  - S0(P0  - Q0)  (19) 


pw(w’v’)o  = - lo 


(20) 


where  F is  the  total  solar  radiation  flux,  B the  infrared 
radiation  flux,  H the  sensible  heat  flux,  LVQ  the  latent 
heat  flux,  P the  precipitation  rate  and  J.  the  wind  stress 
vector.  The  subscript  o indicates  evaluation  of  quantities 
at  the  sea  surface  and  all  fluxes  are  taken  positive  upward. 
Note  that  implicit  in  (20)  is  the  assumption  of  a fully- 
developed  sea. 

The  surface  fluxes  of  sensible  heat,  latent  heat, 
and  momentum  can  be  related  to  mean  field  quantities  using 
the  Monin-Obukhov  theory  (see  Busch,  1977) 

H0  = fJacpc0  (Ur-Uo)(0o-0r)  (21) 


LvQo  = P a^v^q ( ) ( q<}— Qr  ) 

_ _ ^ v 

To  = Pa^m^  ur“Ho ) 


(22) 


(23) 
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where  U is  wind  speed,  X is  wind  velocity,  9 is  potential 
temperature  of  the  air,  and  q is  water  vapor  mixing  ratio 
(mass  of  water  vapor  per  mass  of  dry  air).  The  subscript  r 
indicates  evaluation  of  quantities  at  some  reference  level 
zr  and  the  subscript  o indicates  evaluation  of  quantities 
at  the  sea  surface.  A reference  density  for  air  is  denoted 
by  pa  and  Cp  is  the  specific  heat  at  constant  pressure 
for  air. 
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and  u*  the  surface  friction  velocity  in  the  air  defined  by 
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The  quantity  L in  (25)  is  the  Monin-Obukhov  stability- 
length,  given  to  good  approximation  by 
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Over  most  of  the  ocean,  the  usual  situation  is  one  of  upward 
surface  fluxes  of  sensible  and  latent  boat  which  implies  L<0. 
In  this  case  an  atmospheric  mixed  layer  in  which  S,  q f and  \ 
are  relatively  homogeneous  is  observed  to  exist  extending 
from  the  top  of  the  surface  layer  (at  a height  of  -L)  to 
several  hundred  meters  altitude. 


Eventually  and  in  support  of  OPS,  satellites  will 
routinely  supply  observations  of  mixed  layer  values  for 
potential  temperature,  water  vapor  mixing  ratio,  and  wind 
velocity  over  much  of  the  ocean.  Thus,  it  is  appropriate  in 
this  study  to  relate  the  surface  fluxes  to  mean  field  quanti- 
ties characteristic  of  the  atmospheric  mixed  layer,  which  can 
be  done  by  taking  zr=-L. 


Mon  in-Obukhov  similarity  theory  predicts  that  the 
non-dimensional  gradient  functions  $a  will  be  functions  of 
only  Z/L  in  the  atmospheric  surface  layer.  Proposed  forms 
of  the  <#>  functions  have  been  reviewed  by  Hogstrom  (1974)  and 
it  is  generally  accepted  that  most  surface  layer  data  are 
well  represented  by 


" *q  " - 9 !>’* 


*m  “ - 15  !>”* 


for  L < 0 (31) 
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which  are  the  forms  used  here.  Thus,  making  use  of  the  fact 
that  H0,Q0,  and  tq  are  approximately  constant  in  the 
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atmospheric  surface  layer,  setting  v.r~- L,  and  using  (24), 
(25),  and  (31),  the  bulk  aerodynamic  coefficients  can  be 
written 


CG  = C = 


In 


ft)-  H [inft)  • ° 


.844 


(32) 
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where,  following  many  others  v;e  have  assumed 


zoot  ~ zo  01  = 


(34) 


Furthermore,  following  Clarke  (1970)  we  assume 


z0  = max 


0.032u, 


g 


z , 
min 


(35) 


were  zmin=0.0015  cm.  We  set  $0  equal  to  the  temperature 
at  the  uppermost  gridpoint  of  the  model  ocean  and  thus 
neglect  any  skin  temperature  effects  (see  Grassl,  1976).  We 
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set  qc  equal  to  the  saturation  value  at  the  sea  surface  and 
for  simplicity  we  take  U0  equal  to  zero  and  assume  that  the 

* 

surface  stress  vector  is  in  the  direction  of  the  mixed  layer 
wind  velocity.  Furthermore,  since  we  have  set  zr=-L,  0r,qr, 
and  Ur  are  given  to  excellent  approximation  by  the  values 
of  these  quantities  characteristic  of  the  atmospheric  mixed 
layer  0(„,qm,  and  llm.  Thus,  having  specified  0n,,qra, 
and  V.n  (here  Vm  is  the  wind  velocity  vector  characteristic 
of  the  atmospheric  mixed  layer),  (21) -(23),  (30),  (32),  (33), 
and  (35)  can  be  solved  iteratively  to  yield  H0,  LVQ0,  and  I0. 

Finally  we  calculate  the  surface  radiative  fluxes 
using  the  expressions  presented  by  Wyrtki  (1965) 


F0  = Fl  1 - 0.33(C+C2) 


(36) 


B0  = £b*To4(0.39  - 2.0  q 1 0*)  ( 1 - 0.56C2)  + 


4eb*T03(T0  - T10 ) 


(37) 


where  F^  is  the  surface  solar  flux  in  the  absence  of  clouds, 

C is  the  fractional  cloud  cover,  e the  infrared  emissivity  of 
the  sea  surface,  b*  the  Stef an-Boltzmann  constant,  T0  the  sea 
surface  temperature,  Tio  thc  air  temperature  at  10  m height,  and 
qiO  the  water  vapor  mixing  ratio  at  10  m height.  Since  the 
surface  layer  will  be  thicker  than  10  m for  all  cases  consid- 
ered here  (i.e.,  -L>10  m),  Tjq  and  Qlo  can  '5e  determined  by 
using  (31)  and  integrating  the  profile  relations  (26)  and  (27) 
from  the  surface  to  the  10  m height. 
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DESCRIPTION  OF  THE  EXPERIMENTS 


To  investigate  the  effect  of  inaccuracies  in  the 
meteorological  parameters  (0m,  qm,  Um,  C)  on  a one-dimensional 
forecast  of  upper  ocean  thermal  structure,  a set  of  11  experi- 
ments are  performed.  The  initial  conditions  are  the  same  for 
all  experiments,  characteristic  of  the  summertime  subtropical 
ocean,  and  given  in  Table  1.  The  imposed  meteorological 
parameters  vary  from  experiment  to  experiment,  however,  and 
are  summarized  in  Table  2.  All  model  integrations  are 
carried  out  five  days. 

The  Control  Experiment  (see  Table  2)  can  be  thought 
of  as  a special  case  in  which  the  meteorological  parameters 
are  known  exactly  during  the  course  of  the  integration.  The 
predictions  for  the  other  experiments  will  be  regarded  as 
those  resulting  when  the  meteorological  parameters  are  not 
known  exactly  during  the  forecast  period.  Thus,  anomalies 
between  the  predicted  thermal  structure  for  the  Control 
Experiment  and  the  remaining  experiments  give  measures  of  the 
errors  in  the  ocean  predictions  resulting  from  various 
inaccuracies  ( A 0m , Aqm,  etc.)  in  the  meteorological 
parameters . 

Experiments  1-8  are  designed  to  illustrate  the 
model's  response  to  changes  in  the  meteorological  parameters 
of  standard  magnitude  (i.e.,  1°C,  1 gm  kg-1,  etc.).  Thus,  for 
this  subset  of  the  experiments,  these  quantities  are  varied 
one  at  a time  and  held  constant  during  the  forecast  period. 

In  Experiments  9 and  10  we  simulate  a more  realistic 
forecast  situation  by  allowing  the  inaccuracies  to  occur  in 
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TABLE  1 


Initial  Conditions  for 


all  Experiments 


z(m) 


T(°C) 


S(o/oo ) 


0>z>-20 


27.65 


-20>z>-60  27 . 65+  ( z+20  ) ( 0 . 025  ) 35 


v (ms-1 ) 

-V 


-0. 12.1 


« 


several  meteorolog Leal  parameters  simultaneously  and  grow 
with  time.  The  inaccuracies  increase  during  the  course  of 
the  integration  because  the  errors  in  the  atmospheric  forecast 
will  grow  with  time.  A discussion  of  the  relevance  of  the 
assumed  inaccuracies  for  the  meteorological  paramete.rs  in 
Experiments  9 and  10  can  be  found  in  section  5.2. 

Finally,  we  note  that  values  for  the  empirical  con- 
stants n^-ng  in  the  turbulent  diffusion  model  are  those 
given  by  Mellor  and  Yamada  (1974).  All  other  consants  can  be 
found  either  in  standard  reference  books  or  Table  3.  The 
numerical  scheme  used  is  forward  in  time  with  the  diffusive 
terms  in  ( 1 ) — ( 3 ) treated  implicitly.  A uniform  grid  with  a 
spacing  of  2 m extending  from  the  surface  to  a depth  of  60  m 
is  used.  The  absorption  of  solar  radiation  in  (1)  is  modeled 
using  the  Type  I extinction  profile  from  Jerlov  (1951). 


TABLE  3 


Model  Constants  and  Parameters 


Quintit 


Symbol 


Value 


Coriolis  parameter 


daily-averaged  downward 
flux  of  solar  radiation 
at  sea  surface  for  cloud- 
less conditions 


575  ly  day 


von  Harman' s constant 
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5.1  Model  Response  to  Constant  Inaccuracies 

in  the  Meteorological  Parameters 

In  this  section  we  present  results  from  Experiments  1-8 
in  the  form  of  time-depth  contours  of  temperature  anomaly 
between  the  experiment  in  question  and  the  Control  Experiment. 
These  temperature  anomalies  can  be  regarded  as  errors  in  the 
ocean  forecast  resulting  from  the  constant  prescribed 
inaccuracies  in  the  various  meteorological  parameters 
(A0m,  Aqm,  etc.)  summarized  in  Table  2. 

Fig.  1 shows  the  results  for  Experiment  1 (AOm--l°C). 

The  -2  contour  line  correponds  closely  with  the  base  of  the 
mixed  layer  and  anomalously  cool  water,  caused  primarily  by 
an  increased  flux  of  sensible  heat  from  sea  to  air,  exists 
in  the  mixed  layer.  As  will  be  the  case  for  all  the  experi- 
ments, the  temperature  anomaly  diffuses  rapidly  downward 
from  the  surface  to  within  a few  meters  of  the  mixed  layer 
base,  where  its  downward  diffusion  becomes  much  slower. 

This  is  due  to  the  rapid  decrease  in  turbulent  kinetic 
energy  (and,  consequently,  eddy  diffusivity)  with  depth  near 
the  base  of  the  mixed  layer. 

Note  the  relatively  large  vertical  gradient  in  tempera- 
ture anomaly  at  the  base  of  the  mixed  layer  near  the  end  of 
the  five  day  forecast  (see  Fig.  1).  This  result  will  also 
appear  in  the  other  experiments  and  may  be  significant  since 
the  vertical  temperature  gradient  near  the  mixed  layer  base 
is  an  important  quantity  in  the  theory  of  underwater  sound 
propagation. 
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Also  of  interest  in  Fig.  1 is  the*  fact  that  the  magni- 
tude of  the  temperature  anomaly  in  the  mixed  layer  increases 
more  rapidly  near  the  beginning  of  the  five  day  forecast 
than  near  the  end.  This  occurs,  primarily,  because  the 
mixed  layer  deepens  during  the  period  and,  consequently, 
its  thermal  inertia  is  larger  near  the  end  of  the  experiment. 

Fig.  2 shows  the  results  for  Experiment  2 (Aqm=-1  gm  kg-*). 
In  this  case,  the  anomalously  cool  water  in  the  mixed  layer 
is  due  primarily  to  an  increased  flux  of  latent  heat  from 
sea  to  air.  Comparison  of  Fig.  1 and  Fig.  2 indicates 
remarkable  qualitative  similarity.  The  base  of  the  mixed 
layer  for  Experiment  2 is  again  well  determined  by  the  -2 
contour.  Therefore,  Figs.  1 and  2 show  that  the  time 
evolution  of  the  mixed  layer  depth  is  very  similar  for  both 
Experiments  1 and  2.  Note,  however,  that  the  mixed  layer 
temperature  anomaly  grows  about  10%  faster  for  Experiment  2 
than  for  Experiment  1. 

Fig.  3 shows  the  results  for  Experiment  3 (A  Um=l  ms-*). 
Increased  surface  fluxes  of  both  sensible  and  latent  heat  in 
this  case  are  primarily  responsible  for  producing  anomalously 
cool  water  in  the  mixed  layer.  Note,  however,  that  the 
temperature  anomaly  in  the  mixed  layer  g :ows  slower  in  this 
experiment  than  in  either  Experiments  1 or  2. 


; 


After  day  1,  the  mixed  layer  is  generally  a few  meters 
deeper  in  Experiment  3 than  in  the  previous  two  experiments. 
In  addition,  an  anomalous  quasi-periodic  fluctuation  in 
mixed  layer  depth  with  an  amplitude  of  about  1 m,  which  was 
not  present  in  the  previous  two  experiments,  is  superimposed 
on  the  slow  deepening  rate  of  the  mixed  layer  in  this 


i 


5-3 


I 


experiment.  This  anomalous  fluctuation  in  mixed  layer  depth 
produces  the  quasi-periodic  variation  in  temperature  anomaly 
near  the  base  of  the  layer  evident  in  Fig.  3. 

The  fluctuation  in  mixed  layer  depth  is  associated  with 
the  flux  of  wind  energy  from  air  to  sea  which  is  given  by 
the  dot  product  of  the  surface  wind  stress  vector  and  the 
surface  current  velocity  vector,  T0*jy0.  Although  the 
wind  stress  vector  is  essentially  time  invariant  in  this 
experiment,  the  surface  current  velocity  vector  oscillates 
with  the  inertial  period.  Thus  r0*y0  and,  consequently, 
the  wind  generation  of  turbulent  kinetic  energy  also  varies 
with  the  inertial  period.  Since  the  wind  speed  and,  there- 
fore, the  magnitude  of  the  surface  wind  stress  vector  is 
larger  in  Experiment  3 than  in  the  Control  Experiment,  an 
anomalous  component  is  introduced  in  the  inertial  fluctuation 
of  wind  generation  of  turbulent  kinetic  energy.  This  causes 
the  anomalous  variation  in  mixed  layer  depth. 

Fig.  4 shows  the  results  for  Experiment  4 (AC=0.1).  As 
expected,  the  slightly  increased  cloudiness  produces  a cool 
anomaly  in  the  mixed  layer.  However,  the  magnitude  of  the 
anomaly  produced  is  much  smaller  than  that  resulting  from 
any  of  the  three  previous  experiments. 

Fig.  5 shows  the  results  for  Experiment  5 (Aem=l°C). 

The  warm  anomaly  in  the  mixed  layer  is,  in  this  case,  caused 
primarily  by  a decreased  flux  of  sensible  heat  from  sea  to 
air.  Comparison  of  Fig.  5 to  Fig.  1 shows  that  the  rate  at 
which  the  magnitude  of  the  mixed  layer  temperature  anomaly 
increases  is  faster  in  Experiment  5 than  in  Experiment  1 . 
Thus,  the  response  of  the  mixed  layer  temperature  anomaly  to 


AOm  is  nonlinear.  This  is  caused  partly  by  the  fact  that 
the  mixed  layer  does  not  deepen  as  rapidly  in  Experiment  5 
as  it  does  in  Experiment  1.  For  example,  at  the  end  of  the 
five  day  period,  the  mixed  layer  is  32  m deep  in  Experiment 
1 but  only  26  m deep  in  Experiment  5.  Thus,  more  cold  water 
is  entrained  into  the  mixed  layer  from  below  in  Experiment  1 
than  in  Experiment  5.  In  addition,  the  .average  thermal 
inertia  of  the  mixed  layer  is  greater  in  Experiment  1 than 
in  Experiment  5 since  its  average  depth  is  greater. 

The  dependence  of  the  drag  coefficient  on  static 
stability  of  the  atmospheric  surface  layer  also  contributes 
to  the  nonlinear  response  mentioned  above.  Specifically, 
the  drag  coefficient  for  latent  and  sensible  heat  in  Experi- 
ment 5 is,  on  the  average,  about  15%  smaller  than  that  of 
Experiment  2.  This  models  the  fact  that  the  more  statically 
stable  atmospheric  stratification  in  Experiment  5 tends  to 
reduce  the  intensity  of  the  turbulence  in  the  lower  atmosphere. 

Fig.  6 shows  the  results  for  Experiment  6 
(Aqrn=l  gm  kg“l).  The  warm  anomaly  in  the  mixed  layer 
evident  in  Fig.  6 is  caused  primarily  by  a decreased  flux  of 
latent  heat  from  sea  to  air.  Comparison  of  Fig.  6 to  Fig.  2 
indicates  that  the  response  of  the  mixed  layer  temperature 
anomaly  to  Aqm  is  also  nonlinear.  However,  the  nonlinearity 
is  not  as  large  in  this  case  as  it  was  for  A0m.  This  is 
because  the  drag  coefficient  for  latent  and  sensible  heat  is 
not  as  sensitive  to  A qm  as  it  is  to  A0m  and  is  only  about 
2%  smaller  in  Experiment  6 than  in  Experiment  2.  Thus,  the 
nonlinearity  here  is  essentially  due  only  to  the  difference 
in  deepening  rate  of  the  mixed  layer  between  Experiments  2 
and  6.  The  mixed  layer  depth  in  Experiment  6 behaves  almost 
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exactly  the  same  as  that  of  Experiment  5.  Thus,  at  the  end 
of  the  five  day  forecast,  the  mixed  layer  depth  for  Experiment 
6 is  only  about  26  m while  that  of  Experiment  2 is  32  rn. 

Fig.  7 shows  the  results  for  Experiment  7 
(A  Um=-1  m s-1) . The  warm  anomaly  present  in  the  mixed 
layer  in  this  case  is  due  primarily  to  a decreased  flux  of 
both  sensible  and  latent  heat  from  sea  to  air.  Not  surpris- 
ingly, comparison  of  Kg.  7 to  Fig.  3 indicates  that  the 
response  of  mixed  layer  temperature  anomaly  to  A Um  is  also 
nonlinear.  As  was  the  case  for  the  previous  experiment, 
this  nonlinearity  is  due  almost  entirely  to  differences  in 
mixed  layer  depth  with  differences  in  the  drag  coefficient 
playing  a negligible  role.  Note  that,  after  the  five  day 
forecast,  the  mixed  layer  is  32  m deep  in  Experiment  3 but 
only  26  m deep  in  Experiment  7. 

Also  of  interest  in  Fig.  7 is  the  quasi-periodic  varia- 
tion in  temperature  anomaly  near  the  base  of  the  mixed  layer. 
This  is  a result  of  the  same  mechanism  discussed  previously 
in  connection  with  the  quasi-periodic  fluctuation  in  tempera- 
ture anomaly  evident  in  Fig.  3. 

Fig.  8 shows  the  results  for  Experiment  8 (AC=-0.1). 

The  warm  anomaly  in  the  mixed  layer  is  due  to  the  increased 
downward  flux  of  solar  radiation  at  the  surface.  However, 
as  was  the  case  for  Experiment  4,  the  magnitude  of  the 
anomaly  is  small  compared  to  that  of  the  other  experiments. 
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Model  Response  to  Time-Varying  Inacc uracies 
in  t ho  Meteorol  o gical  Parameto r s 
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In  this  section  we  will  investigate  the  model *s 
response  when  the  inaccuracies  in  the  meteorological  para- 
meters vary  in  a manner  that  simulates  an  actual  operational 
forecast  situation. 

At  the  beginning  of  the  forecast  period,  we 

assume  that  the  magnitude  of  the  inaccuracies  in  the  potential 

temperature,  humidity,  and  wind  speed  of  the  atmospheric 

mixed  layer  are  given  by  nonzero  quantities  5 , <5n>  6 T7. 

b M u 

These  initial  inaccuracies  will  be  present  because  observa- 
tions of  the  initial  atmospheric  state,  which  will  be  made 
by  satellites  over  a large  part  of  the  ocean,  will  not  be 
perfect . 


Estimates  forfig  ,6  q>  and  6 ^ can  be  determined 
from  the  work  of  Moyer  et^  al^.  (1978).  They  compared  data 
deduced  from  the  Nimbus  6 High  Resolution  Infrared  Radiation 
Sounder  and  Scanning  Microwave  Spectrometer  to  data  obtained 
from  weighted  radiosonde  observations  from  the  east  central 
United  States.  They  found  the  standard  deviation  between 
the  satellite  and  radiosonde  observations  of  temperature, 
humidity,  and  geostropliic  wind  speed  at  850  mb,oT»  c?q, 
and  Oy,  to  be  about 

aT  = 1.5  °C 

o - 1.6  gm  kg-1  (38) 

cu  “ 5 m S_1 


» 
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We  will  take 


8 °T 

I °q  <39> 

16  °U 


Thus  6q  and  <S  q are  set  equal  to  one  eighth  standard 
deviation  limits.  Note  that  because  of  surface  friction, 

Um  is  typically  only  about  one  half  the  geostrophic  wind 
speed  in  the  lower  atmosphere.  As  a result,  an  additional 
factor  of  1/2  arises  when  relating  6 g to  a p . Therefore, 
the  value  assigned  to  can  also  be  regarded  as  a one 
eighth  standard  deviation  limit. 

Moyer  et  aT.  (1978)  also  calculated  the  standard 
deviation  between  the  satellite  and  radiosonde  observations 
of  500  mb  geopotential  height,  0$  , and  found 

0$  = 15  m (40) 

This  combined  with  the  results  of  Baumhefner  and  Downey 
(1978)  can  be  used  to  determine  rough  estimates  for  the  rate 
at  which  the  inaccuracies  grow  with  time.  Baumhefner  and 
Downey  did  a forecast  intercomparison  study  for  the  Northern 
Hemisphere  using  six  different  initial  states  and  three 
different  general  circulation  models  of  the  atmosphere.  Data 
from  the  National  Meteorological  Center  operational  analyses 
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were  used  for  forecast  verification.  At  the  end  of  a five 
day  forecast,  the  area-averaged  standard  deviation  between 
the  predicted  and  observed  500  nib  geopotential  height  was 
typically  110  m.  Furthermore,  a good  first  order  approxima- 
tion to  the  results  of  Baumhefner  and  Downey  is  to  assume 
that  the  forecast  error  for  geopotential  height  increases 
linearly  with  time.  Thus,  taking  the  one  eighth  standard 
deviation  limit  for  6^  , we  have 


A$  = 6,  (1  + Yt) 


(41) 


where  A<i>  is  the  inaccuracy  in  500  mb  height,  = 1 and 


110  m 

Y = 15  m - 1 - 1.27  days 

5 days 


-1 


(42) 


Now,  temperature,  humidity,  wind  speed,  and  geopotential 
height  are  expected  to  be  well  correlated  with  one  another. 
Thus,  it  is  reasonable  to  expect  that  the  magnitudes  of  the 
inaccuracies  resulting  from  imperfect  forecasts  of  these 
quantities  will  also  be  well  correlated  with  one  another. 
Therefore,  as  a first  approximation  we  take 


A9  Aq  AU  A$ 
m _ in  _ m _ 


(43) 


;U 


which  implies 

AO  =±60  (1  + Yt) 

£q  = ±<5q  (1  + Yt) 
AU  “tSy  (1  + Yt) 


(44) 


Using  (38),  (39),  (42),  and  (44),  we  find  the  func- 
tional form  for  Aom,  AQm»  and  A^m  given  in  Table  2. 

Note  that  AG  m and  Aqm  are  negative  and  AUm  is  positive 
in  Experiment  9 while  the  reverse  is  the  case  for  Experiment 
10  (see  Table  2).  Thus,  the  inaccuracies  in  Experiment  9 
all  combine  in  such  a way  as  to  increase  the  flux  of  sensible 
plus  latent  heat  from  sea  to  air  while  the  reverse  is  true 
for  Experiment  10. 

The  probability  that  the  inaccuracies  in  the  atmos- 
pheric parameters  for  a real  forecast  situation  will  have 
the  same  signs  as  and  magnitudes  as  least  as  large  as  those 
prescribed  for  Experiments  9 and  10  can  be  estimated  as  fol- 
lows. First,  we  define  the  inaccuracies  in  temperature, 
humidity,  and  wind  speed  for  a hypothetical  forecast  situation 
to  be  AG*,  A q*,  and  A U* , respectively,  and  assume  that 
they  are  normally  distributed  about  zero.  In  addition  we 
assume  that  the  standard  deviations  of  these  quantities  at 
the  beginning  of  the  forecast  period  on  some  Northern  Hemi- 
sphere grid  are  given  by  our  earlier  chosen  values,  , 6 q, 
and6y  Now,  if  the  standard  deviations  for  AG  *,  Aq*, 
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and  A U*  increase  at  least  as  fast  with  time  during  the  fore- 
cast period  as  the  standard  deviation  of  AO  based  on  the 
linear  fit  to  the  results  of  Baurnhefner  and  Downey  (1978), 
then  our  prescribed  values  for  AG  m,  A qni,  and  Um  (see 
eq.  44)  will  always  lie  within  the  one  eighth  standard 
deviation  limits  for  AO  *,  A q*,  and  A U* . 

Now,  to  obtain  an  estimate  of  the  lower  bound  on 
the  probability  in  question,  we  will  assume  for  the  moment 
that  AO  *,  A q*,  and  U*  are  completely  uncorrelated 
with  one  another.  From  probability  theory  this  implies  that 
AO  *,  A q*,  and  A U*  are  statistically  independent  which 
means  that  the  joint  probability  distribution  function  for 
these  quantities  can  be  written  as  a product  of  the 
individual  probability  distribution  functions.  Therefore, 
based  on  the  one  eighth  standard  deviation  limits,  the 
probability  that  the  magnitudes  of  AO  *,  A q*,  and  A u* 
will  be  simultaneously  larger  than  the  magnitudes  of  our 
prescribed  values  for  AO  m,  A qm,  and  A Um  is  93%. 

However,  the  probability  that  AO*,  A q* , and  A U*  will 
also  have  the  same  signs  as  the  prescribed  values  in 
Experiment  9 is  only  93%/8=12%,  as  is  also  the  case  for 
Experiment  10.  Thus,  even  if  the  inaccuracies  in  the 
atmospheric  parameters  are  completely  uncorrelated,  we 
expect  that  they  will  combine  to  yield  inaccuracies  in  the 
surface  flux  as  least  as  large  as  those  in  either  Experiments 
9 or  10  over  about  24%  of  the  Northern  Hemisphere  ocean. 

In  reality,  however,  we  expect  that  the  inaccuracies 
in  the  atmospheric  parameters  will  be  correlated  with  one 
another.  For  example,  AO*  and  A q*  should  be  very  well 
correlated  with  a positive  correlation  coefficient  as 


demonstrated  by  the  following  qualitative  argument.  Both 
temperature  and  humidity  in  the  lower  atmosphere  generally 
decrease  from  equator  to  pole.  Therefore,  if  at  some 
location  during  the  course  of  the  forecast  period  the  real 
wind  in  the  lower  atmosphere  develops  a stronger  poleward 
component  than  the  forecast  wind,  horizontal  advection  will 
tend  to  force  both  AG*  and  Aq*  negative.  The  reverse 
case  will  occur  if  the  real  wind  develops  a stronger  equator- 
ward  component  than  the  forecast  wind,  and  AG  * and  A q* 
will  again  be  positively  correlated.  Correlation  between 
inaccuracies  in  the  temperature  and  inaccuracies  in  the  wind 
speed  is  less  obvious.  This  is  because  the  inaccuracies  in 
the  wind  field  are  related  to  horizontal  gradients  of  the 
inaccuracies  in  the  temperature  field,  rather  than  the 
temperature  inaccuracies  themselves. 

Therefore,  as  a rough  estimate  of  the  upper  bound 
on  the  probability  that  the  inaccuracies  will  be  at  least  as 
bad  (from  an  ocean  forecasting  point  of  view')  as  those  pre- 
scribed in  either  Experiments  9 or  10,  v/e  here  assume  that 
AG  * and  A q*  are  perfectly  correlated  with  each  other  but 
completely  uncorrelated  with  AD*.  From  arguments  similar 
to  those  presented  earlier  in  this  section,  we  find  an 
estimate  of  this  upper  bound  to  be  45%.  Thus,  to  summarize 
this  analysis,  we  expect  the  errors  in  real  ocean  forecasts 
due  to  inaccuracies  in  the  air  temperature,  humidity,  and 
wind  speed  alone  to  be  at  least  as  bad  as  either  those 
resulting  from  Experiments  9 and  10  over  roughly  24%  to  45% 
of  the  ocean  in  the  Northern  Hemisphere. 

Fig.  9 shows  the  results  for  Experiment  9.  In  this 
case,  the  mixed  layer  becomes  anomalously  deep  and  cool 
because  of  increased  wind  generation  of  turbulent  kinetic 
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energy  in  the  upper  ocean  and  increased  fluxes  of  sensible 
and  latent  heat  from  sea  to  air.  Note  that  the  mixed  layer 
temperature  anomaly  grows  more  rapidly  near  the  end  of  the 
five  day  forecast  tha*1  near  the  beginning.  This  is  because 
the  atmospheric  prediction  becomes  progressively  worse  with 
time  during  the  forecast  period  (and  hence  the  magnitudes  of 
A6m,  Aq,n,  etc.  become  progressively  larger). 

At  the  end  of  the  five  day  period,  the  mixed  layer 
is  8 m deeper  than  that  of  the  Control  Experiment . This 
accounts  for  the  warm  anomaly  near  the  base  of  the  mixed 
layer  shown  in  Fig.  9 and  represents  a significant  forecast 
inaccuracy.  Also  of  interest  on  Fig.  9 is  the  region  of 
large  vertical  gradient  of  temperature  anomaly  between  28 
and  38  m depth  at  day  five.  As  mentioned  earlier,  this 
concentration  of  the  isolines  of  temperature  anomaly  near 
the  base  of  the  mixed  layer  may  have  important  implications 
for  acoustic  prediction. 

Fig.  10  shows  the  results  for  Experiment  10.  The 
mixed  layer  becomes  anomalously  shallow  and  warm  in  this 
experiment  because  the  wind  generation  of  turbulent  kinetic 
energy  in  the  upper  ocean  and  the  surface  fluxes  of  latent 
and  sensible  heat  are  anomalously  small.  As  was  the  case 
for  Experiment  9,  the  mixed  layer  temperature  anomaly 
increases  faster  near  the  end  of  the  forecast  than  near  the 
beginning.  Furthermore,  since  the  average  thermal  inertia 
of  the  mixed  layer  is  anomalously  small  in  this  case  (because 
of  its  shallowness)  the  effect  is  even  more  pronounced  than 
it  was  in  Experiment  9. 

Tn  addition,  comparison  of  Figs.  9 and  10  shows 
that  the  mixed  layer  temperature  anomaly  resulting  at  the 
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end  of  the  five  day  forecast  is  about  twice  as  large  in 
Experiment  10  than  in  Experiment  9.  Since  the  inaccuracies 
in  the  atmospheric  parameters  (A0,n,  Aqra,  and  A Ujn ) 
differ  between  Experiments  9 and  10  in  sign  but  not  in 
magnitude  (see  Table  2),  this  reflects  the  nonlinci.rities 
discussed  in  section  5.1. 

At  the  end  of  the  five  day  forecast,  the  predicted 
mixed  Layer  depth  for  this  experiment  is  18  in  shallower  than 
that  of  the  Control  Experiment.  This  is  a severe  forecast- 
inaccuracy  and  accounts  for  the  region  of  large  vertical, 
gradient  in  temperature  anomaly  between  10  and  30  m depth. 
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6.0  SUMMARY  AND  CONCLUSIONS 

We  have  used  a one-dimensional  model  of  the  upper 
ocean  to  investigate  the  effect  of  inaccuracies  in  the 
potential  temperature,  specific  humidity,  and  wind  speed  of 
the  atmospheric  mixed  layer  and  the  fractional  cloud  cover 
(which  are  all  regarded  as  imposed  meteorological  parameters) 
on  forecasts  of  ocean  thermal  structure.  Vertical  eddy 
fluxes  of  heat,  salinity,  and  momentum  in  the  oceanic  mixed 
layer  are  parameterized  using  the  Level-2  turbulence  closure 
model  of  Mel lor  and  Yamada  (1974).  Surface  eddy  fluxes  of 
sensible  heat,  latent  heat,  and  momentum  are  parameterized 
with  bulk  aerodynamic  coefficients  determined  from  Monin- 
Obukhov  similarity  theory.  Surface  fluxes  of  solar  and 
infrared  radiation  pro  parameterized  with  the  formulas  of 
Wyrtki  (1965). 


A Control  Experiment  is  performed  by  imposing  time- 
invariant  values  for  the  meteorological  parameters  and 
integrating  the  model  ocean  forward  in  time  for  five  days 
from  an  initial  state  characteristic  of  the  summertime 
subtropical  ocean.  A set  of  ten  additional  experiments  are 
performed  by  integrating  the  model  for  five  days  from  the 
same  initial  state  used  for  the  Control  Experiment  but  with 
different  values  for  the  various  imposed  meteorological 
paramet  ers . 

The  predicted  upper  ocean  thermal  structure  for  the 
Control  Experiment  is  regarded  as  a forecast  in  which  the 
meteorological  parameters  are  known  exactly  during  the  fore- 
cast period.  The  predicted  thermal  structures  for  the  remain- 
ing ten  experiments  are  considered  to  be  those  resulting  when 
these  parameters  are  not  known  exactly  during  this  period. 


s 


Time-depth  contours  of  temperature  anomaly  between  each  of 
the  ten  experiments  and  the  Control  Experiment,  which  can  be 
thought  of  as  the  error  introduced  into  the  ocean  forecast 
by  inaccuracies  in  the  meteorological  parameters,  are 
presented  and  discussed. 


In  all  cases,  the  temperature  anomaly  diffuses 
rapidly  downward  from  the  surface  to  within  a few  meters  of 
the  mixed  layer-  base.  Because  of  the  marked  decrease  in 
eddy  diffusivity  there,  downward  diffusion  of  the  anomaly  is 
much  slower  below  this  level.  This  results  in  a relatively 
large  vertical  gradient  in  temperature  anomaly  near  the  base 
of  the  mixed  layer  which  may  have  important  implications  for 
acoustic  modeling.  Both  the  depth  of  the  mixed  layer  and 
the  magnitude  of  the  mixed  layer  temperature  anomaly  are 
found  to  be  nonlinear  functions  of  the  inaccuracies  in  the 
potential  temperature,  humidity,  and  wind  speed.  The 
magnitude  of  the  anomalies  in  mixed  layer  depth  and  mixed 
layer  temperature  resulting  from  an  inaccuracy  in  fractional 
cloud  cover  of  0.1  is  small  compared  to  the  anomalies 
produced  by  inaccuracies  in  the  temperature,  humidity,  and 
wind  speed  of  1°,  1 gin  kg"1,  and  1 m s_1 , respectively. 


Estimates  of  the  time  evolution  of  the 
inaccuracies  in  the  air  temperature,  humidity,  and  wind  speed 
during  an  actual  operational  forecast  situation  are  made. 
Expected  magnitudes  of  the  inaccuracies  at  the  beginning  of 
the  forecast  period  are  obtained  from  the  satellite/radiosonde 
comparison  study  of  Moyer  et_  al_.  (1978).  This,  combined 
with  the  numerical  weather  prediction  study  of  Baumhefner 
and  Downey  (1978),  is  used  to  estimate  the  rate  at  which  the 
magnitude  of  the  inaccuracies  grow  with  time  during  the  fore- 
cast period.  Then,  when  the  sign  of  the  various  inaccuracies 
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t are  chosen,  rough  estimates  of  the  probability  that  the 

prescribed  inaccuracies  will  be  characteristic  of  an  actual 
forecast  situation  are  made.  In  both  Experiments  0 and  10, 
for  which  the  inaccuracies  are  prescribed  in  this  manner, 

* this  probability  estimate  turns  out  to  bo  about  12-22%. 

In  Experiment  9 during  the  entire  forecast  period, 
the  signs  of  the  inaccuracies  are  such  that  the  air  temperature 
and  humidity  are  anomalously  low  and  the  wind  speed  is  anoma- 
lously high.  At  the  end  of  a five  day  forecast  a mixed  layer 
depth  anomaly  of  8 m and  a mixed  layer  temperature  anomaly 
of  -0.2°C  results.  In  Experiment  10,  the  inaccuracies  in 
the  three  atmospheric  parameters  each  have  the  same  magnitude 
as  those  in  Experiment  9 but  opposite  sign.  In  this  case, 
the  mixed  layer  depth  anomaly  and  temperature  anomaly 
resulting  from  a five  day  forecast  are  -18  m and  0.38°C 

* respectively. 

Of  course  these  results  are  not  general  since  the 
magnitude  of  the  anomalies  in  mixed  layer  depth  and  tempera- 

* ture  produced  by  the  inaccuracies  in  the  forecast  atmospheric 
parameters  depend  on  the  initial  depth  of  the  mixed  layer 
and  the  stratification  below  it.  For  example,  during  the 
winter  when  the  mixed  layer  may  be  over  100  m deep,  the 

* magnitude  of  the  anomalies  may  be  only  1/3  to  1/5  of  those 
found  here.  Nevertheless,  the  minimum  prediction  performance 
criteria  for  OPS  requires  an  accuracy  in  temperature  of 
0.5°C  (Grabowski  and  Roberts,  1978).  Thus,  the  present 

? study  indicates  that,  at  least  when  the  mixed  layer  is  shallow, 

errors  introduced  by  inaccuracies  in  the  meteorological  forc- 
ing  alone  can  amount  to  a substantial  part  of  the  minimum 

prediction  performance  criterion. 

l 
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Errors  introduced  into  ocean  forecasts  due  to 


inaccuracies  in  the  atmospheric  forcing  will  be  largest  in 
the  mixed  layer.  Fortunately,  because  of  satellite  observa- 
tions of  sea  surface  temperature,  the  initial  inaccuracies 
in  mixed  Layer  temperature  for  an  ocean  forecast  will  be 
relatively  small.  Thus,  a trade  off  Situation  exists. 
Inaccuracies  in  the  atmospheric  forcing  will  introduce 
errors  into  forecasts  of  the  upper  ocean,  but  the  initial 
conditions  for  the  forecast  will  be  relatively  good  there. 
Inaccuracies  in  the  atmospheric  forcing  will  have  little 
effect  on  the  forecast  state  of  the  deep  ocean,  but.  the 
initial  conditions  for  the  forecast  will  be  relatively  bad 
there . 


The  present  study  can  be  extended  in  several  ways. 
Numerical  weather  prediction  studies  using  data  routinely 
generated  by  an  operational  forecast  model  and  analysis 
program  could  be  done  with  particular  regard  for  the  fields 
that  are  important  to  the  ocean  forecast  problem.  In 
addition,  spatial  and  seasonal  variability  of  the  atmospheric 
forecast  errors  might  be  obtained.  For  example,  the  errors 
may  be  generally  small  in  the  trade  wind  region  compared  to 
those  of  the  midlatitudes,  where  they  may  be  especially 
large  in  winter. 

Upwelling  and  downwelling  are  induced  in  the  mid 
ocean  by  surface  wind  stress  curl.  Thus,  forecast  errors  of 
wind  stress  curl  should  also  be  analyzed.  If  these  errors 
prove  to  be  large  enough,  then  a study  similar  to  the 
present  one  should  be  performed  with  a three-dimensional 
model . 
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APPENDIX 
LIST  OF  SYMBOLS 

a thermal  expansion  coefficient  for  seawater 

Aj{,  Ajj  stability  functions  Tor  eddy  coefficients  in 
oceanic  mixed  layer 

b saline  expansion  coefficient  for  seawater 

b*  Stefan-Boltzmann  constant 

B()  surface  flux  of  infrared  radiation 

c,  Cp  specific  heat  at  constant  pressure  for  seawater, 
air 

C fractional  cloud  coverage 

drag  coefficient  for  sensible  heat,  latent  heat, 
momentum 

E turbulent  kinetic  energy 

f Coriolis  parameter 

F downward  flux  of  solar  radiation 

F0>  Fi  downward  flux  of  solar  radiation  at  sea  surface, 
daily-averaged  downward  flux  of  solar  radiation 
at  sea  surface  for  cloudless  conditions 

g acceleration  of  gravity 

II0  • surface  flux  of  sensible  heat 

Kjj,Kg,K..  eddy  diffusion  coefficients  for  heat,  salinity, 
and  momentum 

I.  Monin-Obukhov  stability  length 

Ly  latent  heat  of  evaporation  for  water 

1 turbulence  length  scale 

nl-n8  empirical  constants  in  eddy  diffusion  model 
P0  precipitation  rate  at  sea  surface 

q,  qni  water  vapor  mixing  ratio,  water  vapor  mixing 

ratio  of  the  atmospheric  mixed  layer 

Qo^lr’^lO  wa^er  vapor  mixing  ratio  at  sea  surface,  reference 
height,  ten  meters  height 

surface  evaporation  rate 


flux  Richardson  number 
sal  ini  ty 

reference  salinity,  salinity  at  sea  surface 
temperature 

reference  temperature,  temperature  at  sea  surface 
temperature  at  ten  meters  height 


time 

wind  speed  of  atmospheric  mixed  layer 

wind  speed  at  sea  surface,  wind  speed  at  reference 
height 

surface  friction  velocity  in  air 

wind  velocity  vector  at  reference  height 

wind  velocity  vector  of  the  atmospheric  mixed 
layer 

current  velocity  vector,  cut  ;nt  velocity  vector 
at  sea  surface 

Cartessian  coordinates,  x positive  eastward,  y 
positive  northward,  z positive  upward  from  the 
sea  surface 

unit  vectors  in  the  x,y,z  directions 

reference  height,  surface  roughness  length 

infrared  emissivity  of  sea  surface 

initial  inaccuracies  of  potential  temperature, 
water  vapor  mixing  ratio,  wind  speed,  geopotential 
height 

von  Harman's  constant 

standard  deviation  of  inaccuracies  in  temperature, 
water  vapor  mixing  ratio,  wind  speed,  geopotential 
height 

density 

reference  density  for  water,  air 

potential  temperature,  potential  temperature  of 
atmospheric  mixed  layer 

potential  temperature  at  sea  surface,  potential 
temperature  at  reference  height 


t 


<**  geopotential  height  of  500  mb  surface 

4>q.  <?>  , 4>m  nondimensional  vertical  derivatives  of 

1 1 temperature,  water  vapor  mixing  ratio, 

and  wind  speed  in  the  atmospheric  sur- 
face layer 


-To 

surface  wind  stress  vector 

AB  , 
m 

An  AU  , AC 
in’  m 

prescribed  inaccuracies 
and  C 

in  o ,. 
m ’ 

q , U , 
‘m  m’ 

* 

AG  , 
m’ 

AV  Aum 

inaccuracies  in  0 , q . 
hypothet i cal  forecast 

and  U 

m 

for  a 

A$ 

inaccuracy  in  <I> 

(“)  spatial  average  at  constant  z taken 

across  a region  large  enough  to  encom- 
pass a statistically  significant  sam- 
pling of  all  unresolved  r . nomena 


( ' ) departure  from  above  defined  spatial 

average 
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